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ABSTRACT 

We present an analysis of the first space-based far-IR-submm observations of M 33, which measure the emission from the cool 
dust and resolve the giant molecular cloud complexes. With roughly half-solar abundances, M33 is a first step towards young 
low-metallicity galaxies where the submm may be able to provide an alternative to CO mapping to measure their H2 content. In 
this Letter, we measure the dust emission cross-section a using SPIRE and recent CO and Hi observations; a variation in g is 
present from a near-solar neighborhood cross-section to about half-solar with the maximum being south of the nucleus. Calculating 
the total H column density from the measured dust temperature and cross-section, and then subtracting the H I column, yields a 
morphology similar to that observed in CO. The H2/H1 mass ratio decreases from about unity to well below 10% and is about 
15% averaged over the optical disk. The single most important observation to reduce the potentially large systematic errors is to 
complete the CO mapping of M 33. 
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1. Introduction, data, and dust temperature 

Understanding star formation requires studying the inter- 
play between the phases of the interstellar medium (ISM) . 
Dust processes most of the energy transiting the ISM, but 
the cool dust component, although representing the vast 
majority of the dust mass, is difficult to observe from 
the groimd. Herschel SPIRE observations (jPilbrattl [20Tot 
I Grif fin 2010) arc the first space-based 250-500 ^m data and 
as such provide a unique occasion to put together a global 
picture of the cool gas and dust in M 33. In particular, we 
compare the morphology of the Far-IR emission and that 
of the gas as determined from CO and Hi measurements 
and attempt to measure how the dust cross-section varies 
in M 33. A longer term goal is to be able to use the dust 
emission to constrain the variation of the A^(Il2)//co factor 
within M 33 and elsewhere. 

This Letter is one of a series on the HERM33ES project 
on the ISM of the Local Group galaxy M 33, an overview of 
which is given in iKramer et al.l (|2010[ 1. hereafter KIO. For 
consistency with the other M 33 papers in this volume, we 
adopt a distance of -D = 840 kpc for M33 (i.e., 25 arcsec = 
100 pc) and orientation parameters of PA = 22.5 degrees 
and i = 56°. Wc use the recent HERM33ES SPIRE observa- 
tions at 250, 350, and 500 /im combined with C0(2~l) ob- 
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servations from the M33CO@IRAM project, as a tracer of 
the molecular component , and a high- resolut ion mosaic of 
VLA HI data (both from lGratier et al.ll2010[) . The SPIRE 
data were first processed as described in KIO and then con- 
verted from Jy/beam to brightness units (MJy/sr). To esti- 
mate the dust temperature, the 250/xm data were convolved 
to the 350/im beamsize (~ 25") and, assuming a single tem- 
perature grey body with an emissivity oc v^ with /3 = 2, a 
temperature was derived from the flux ratio. At the tem- 
peratu res of the cool component seen in M 33. lDuDac et al.l 
(I2003D find /3 w 2. To minimize the effect of the uncer- 
tainty in /3, we chose to use adjacent bands to estimate 
temperature. The 250/350//m ratio provides more accurate 
temperatures than the 350/500 micron ratio even for tem- 
peratures below 10 K. A 15% variation (or uncertainty) in 
the 250/350/im ratios corresponds to a temperature change 
of 0.9, 2.1, and 4.1 K at temperatures of 10, 15, and 20 K 
but the same uncertainty in the 350/500/im ratios yields 
1.3, 3.2, and 6.3 K errors for the same dust temperatures. 
A further advantage is that we obtain the dust temperature 
at a resolution typical of giant molecular clouds (GMC), 
- 100 pc. 

Ideally, a multi-component fit would be used but this re- 
quires high S/N data at many wavelengths over the whole 
disk. By assuming that the cool dust component domi- 
nates the emission beyond 250/zm, we could calculate the 
dust temperature out to the optical radius of M 33. The 
resulting temperature map is shown in Fig. 1. Figure 2 
shows a comparison with the temperature of the cool com- 
ponent of the preferred two-component (warm plus cool 
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Fig. 1. Dust color temperature map of M 33 calculated 
from the SPIRE 250 and 350 ^\n flux ratio with /3 = 2 at 
the ~ 25" resolution of the 350/im data. Coordinates are 
J2000, the temperature scale is shown at the top in Kelvin, 
and the ellipse shows the 8 kpc radius, slightly beyond R25. 



dust) model fit to data between 24^m and 500//m from 
KIO. The two-component model uses /? = 1.5 and the tem- 
peratures are higher out to 6 kpc. The temperatures in 
the radial bins have been estimated by averaging the tem- 
peratures rather than averaging the emission as in KIO. 
The latter yields slightly higher temperatures because the 
dust is usually warmer where emission is strong (compare 
Fig. 1 with Fig 1. in KIO). KIO calculate the dust col- 
umn density using k — 0.4(i'/250GHz)^ cm^g"-'^ of dust, 
which is equivalent to a dust cross-section per H-atom of 
cr = 1.1 X 10^^^(A/250/im)~^cm^ for a hydrogen gas-to- 
dust mass ratio of 140. In the following we use the SPIRE 
250/350 /im color temperature because the data cover a 
greater area and agree well with the cool dust temperature 
of the two-component fit, showing the domination of the 
cool dust at these wavelengths. No correction for line con- 
tamination was subtracted from the SPIRE data. At these 
frequencies, the CO lines contribute very littl e to the con- 
tinuu m fiux, unlike at 1.3mm or 850/.tm fe.g. iBraine et aLl 
[T99I . 

Figure 3 shows the SPIRE 250/im emission with the 
C0(2-l) emission as contours at 25" resolution to show how 
closely the CO emi ssion follows the dust emission peaks. 
iGratier et al.l (|2010( ) show a detailed comparison of the CO 



and H I emission with star formation tracers such as Ha 
and Spitzer 8 and 24/im maps. The 250/im-bright regions 
are detected in CO and the general morphology of the cool 
dust emission is seen in the Hi image. 
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Fig. 2. Radial distribution of the average dust tempera- 
ture in the stellar disk of M33. The solid line joining the 
triangles shows the 250/350 yum color temperature and for 
comparison we show the cool dust temperature (large dots) 
obtained by KIO in their preferred two-component model 
with /3 = 1.5. The error bars on the triangles indicate a 15% 
calibration uncertainty in the 250/350 fira ratio converted 
to temperature. 



2. Dust cross-section 

In order to estimate the total gas mass of a galaxy from 
its dust emission, it is first necessary to measure the dust 
emission cross-section per H-atom a. We do this by se- 
lecting regions with HI emission and well-constrained dust 
temperatures but with little or no CO emission (in order 
to avoid the assumption of a N(H2)/Ico value and CO 
line ratio). Thereby, we can relate the observed neutral 
hydrogen column density, Nhi, to the dust optical depth, 
T = (tNh = S^/B^^T, and infer the cross section a. In 
practice, two methods were used: (i) we selected large con- 
tiguous regions as shown on the map in Fig. 3 and (m) all 
pixels with HI emission, a defined dust temperature, and 
no detected CO emission (Ico < 0.05 K km s~^) were taken, 
using maps at the resolution of the 350/im SPIRE data. In 
the optically thin limit. 



A^tot = iVffl + Ah, = S,/ia B,,t) 



(1) 



where Sn is the flux density and A'tot is the H column den- 
sity. Where CO is not detected (i.e. A'tot ~ A'hi), we can 
calculate a from the dust temperature and the observed H i 
column density and SPIRE flux by cr w 5'i//(A^hi B^t)- 

From the regions indicated in Fig. 3, we estimate a (us- 
ing Eq. 1) for different radii and the results are shown in 
the flrst line of Table 1. We assume that the cross section 
within the inner 4 kpc is the same as at 4kpc for two rea- 
sons: {i) there are few truly CO-free r egion s in the central 
part of M 33 and (ii) iGardan et al] (|2007t ) showed (their 
Fig. 13) that the H1-H2 relation was similar within 4kpc 
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Fig. 3. SPIRE 250/^m emission with the IRAM C0(2-l) 
emission at 25" resolution superposed as contours and with 
the polygons to the upper left used to calculate the dust 
cross-section per H-atom. The beams are shown as white 
dots in the lower left corner. The CO contours are at 0.5 
(red), 2 (blue), and 4 K km s"^ in the main beam scale. 
The large polygon indicates the area covered so far by the 
IRAM CO survey. 



but beyond 4 kpc the H2 fraction decreased sharply even 
at constant total Hydrogen column density, so we consider 
R ^ 4kpc as a transition in ISM properties. It is very likely 
that some molecular gas is present in these regions, par- 
ticularly those closer to the center, so a is likely overesti- 
mated. In the second approach, in order to exclude averag- 
ing in undetected molecular clouds, we take the peak of the 
histogram of a values for the areas without detected CO 
emission but with HI emission and a constrained dust tem- 
perature, yielding the lower values for cr, as shown in Table 
1. Using this method, we find a north-south difference in 
M33 at equivalent radii. The uncertainty in the histogram 
method is about ±0.1; although the distribution is some- 
times broad, the peak value is well defined. The histograms 
of a values are clearly different in the north and south. The 
polygon method was difficult to use in the south due to the 
absence of large CO-free regions with HI and dust emission. 
An advantage of the histogram method is that it excludes 
the tail of high a which may be due to H2 without detected 
CO emission. 



The intrinsic expectation is that because the oxygen 
abundance i s about half solar, t he dust cross-section should 
be as well (jDraine et al.l [20071) . a nd that given the shal- 
low abundance gra dient in M33 (jRosolowskv et al] 120071 : 
^ Magrini et al.ll2009l ). this should hold for the entire galaxy. 
However, it quickly became apparent that the dust cross- 
section per H-atom varied over the galaxy, decreasing with 
radius by close to a factor 2, with a higher abundance in 
the so uth - similar to the variation found by iMagrini et al.l 
(l2010f ) from optical HII region and PNe abundance mea- 
surements. In the solar n eighborhood, a is about a k, 1.1 x 
lO^^cm^ per H-atom (jDraine fc Led Il984t iDraine fc LI 
l 2007[ ) at 250^m and varying as v^ . If we apply a cr of half 
this value to M33 using the dust temperature map in Fig. 1, 
we obtain unrealistically high gas masses in the inner disk: 
^ 40 M0 pc~^ on average, wherea s the H i is about 1 1 M© 
pc^^ and the H2 considerably less (jGratier et al.ll2010l) . For 
a half-solar a, the (cool) dust temperature would have to 
be about 30K, beyond any experimental uncertainties, so 
we conclude that a must be higher in the inner disk. In 
the outer disk, however, the H i column density (where no 
CO emission is present) is roughly equal to what is derived 
from the dust with a = 0.5(70, suggesting that this value is 
appropriate for the outer disk. 



Table 1. Dust cross-section a at 250/im as a function of 
radius in M33, expressed in units of 1.1 x 10~^^cm^ per 
H-atom. The polygons are shown in Fig. 3. The histogram 
method was applied separately in the north and south (lines 
2 and 3), leading to the "Model" values used to estimate 
the total H column density to make Fig. 4. 
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3. Total and molecular gas mass 

To really obtain an accurate gas mass from the cold dust 
emission data, the calibration uncertainties in the SPIRE 
bands must be reduced because the 15% uncertainty is not 
sufficient to distinguish between different solutions of /3 and 
temperature. A second requirement is that the whole disk 
of M33 must be observed in CO so that a well-determined 
variation of a^ based on the second technique, can be con- 
structed. The Hi data provide a reliable picture of the Hi 
column density, so subtracting this from the total gas col- 
umn density derived from the dust emission yields an es- 
timate of the molecular gas mass which does not rely on 
CO emission. A further caveat is that a must not change 
from the atomic to molecular medium. Once the [C i] and 
[C 11] lines have been measured, these uncertainties can be 
reduced by the additional constraint that the C/H ratio 
(from the sum of the C reservoirs) should be equal to the 
C abundance derived by other means. 

Using the information currently available, we have used 
Table 1 (second method) to build an image of the dust 
cross-section a, erring on the low-cr side, as indicated in the 
last line of Table 1. We then used this a, scaled to 500/im 
with /3 = 2, and the dust temperature (Fig. 1) to estimate 
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Fig. 4. H2 column density map estimated from the 500^m 
SPIRE map, /? = 2, the 250/350At m dust temperature a s 
shown in Fig. 1, the Hi data from iGratier et all (J2010t) . 
and the cross-sections given in Table 1. The column density 
scale is shown at the top in units of H2 molecules per cm^. 
CO contours at 0.4 and 1.6 K km/s (black) and 5 K km/s 
(blue) are superposed. 



the total gas mass from the 500//m data and then sub- 
tracted the H I column density to obtain the H 2 column den- 
sity u sing the interferometric Hi maps from iGratier et aTl 
(2010"), which recove r more than 90% of the flux found by 
[Putman et al.l (|2009l ) using Arecibo. The missing 21cm flux 
is expected to be located in the mid-to-outer disk where the 
rotation curve is fairly flat. A morphological comparison of 
the dust-derived H2 column and the CO contours in Fig. 4 
shows that while major regions were not missed by the CO 
observations, the dust-derived H2 map does not correlate 
perfectly with the CO. For example, the strong CO peak at 
(01^34'"09".4,-h30°49'06") corresponds to a red (not white) 
region in Fig. 4; this may of course reflect a real variation 
in the A^(H2)//co ratio. 

The dat a used in the sca tter plot in Fig. 5 (see also 
Fig. 7 in iLerov et"aLl |2009[ ) yield a formal fit with a 



A^(H2)//C0{2- 



1) 



3 X 10^° H2 per K km s ^ value, corre- 



sponding to a roughly galactic iV(H2)//co(i-o) = 2.1 x 10 " 
H2 per K km s^^ for a line ratio C0(y5h) = 0.7. The inner 
disk A^(H2)//co(2-i) is lower and rises to about 4 x 10^° 
H2 per K km s~^ beyond 2 kpc. The offset, in principle. 



4 r 



E 3 



2 r 



i- \ 


. / ■ ■ . 


\ ]^v 


y 


^'^■^,un.-s- ■: 


- X ^ 


■^^gm-^-- : 


.- ij^ ■' ^^ 




m 


W" - 




\..>M 




: 














: 1 
= . , . 1 . 


-. 



-202468 
C0(2-l) (K km s"') 

Fig. 5. Scatter plot of H2 column density as derived from 
the dust emission and the CO (2-1) intensity. The solid red 
line is a fit to the entire dataset and corresponds to a 
iV(H2)//co(2-i) value of 3 x lO^o H2 per K km s^^ with 
an offset as described in the text. A radial variation is 
present: fitting within 2 kpc yields a A^(H2)//co(2-i) value 
of 2.2 X lO^o H2 per K km s'^ and 4.1 x lO^o H2 per K km 
s~^ beyond 2 kpc (dotted and dashed lines, respectively). 



implies that some low column density H2 is not seen in 
CO. Both values depend strongly on the details of the dust 
cross-section, illustrating the need for a complete map of M 
33 in CO. 

The total hydrogen mass derived from the dust emission 
within 8kpc is estimated to be 1.63x10^ Mq, very similar to 
the m ore directly measured H 1-I-H2 masses in lGratier et al.l 
(|2010[ ). The dust cross-sections we derive suggest dust-to- 
gas mass ratios (including He) ranging from ^ 125 in the 
inner 4 kpc to 200 near R25 for Milky Way like dust (c/. 
KIO). Compari ng the total dust-de rived H2 mass with the 
CO luminosity (jGratier et al.ll2010D yields a A^(H2)//co of 
about 1.5 times the Galactic A^(H2)//co- Because the Hi 
dominates, however, a small error in the total dust-derived 
gas mass translates into a large uncertainty in the dust- 
derived H2 mass. It is clearly necessary to obtain whole- 
galaxy measurements of a in CO-frce zones to refine these 
estimates. 
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